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Abstract: A fast transition between a standard matter- like era and a late ACDM-like 
epoch (or more in general, a CDM+DE era), generated by a single Unified Dark Matter 
component, can provide a new interesting paradigm in the context of general relativity, 
alternative to ACDM itself or other forms of DE or modified gravity theories invoked to 
pH \ explain the observed acceleration of the Universe. UDM models with a fast transition 

have interesting features, leading to measurable predictions, thus they should be clearly 
distinguishable from ACDM (and alternatives) through observations. Here we look at 
different ways of prescribing phenomenological UDM models with fast transition, then 
focusing on a particularly simple model. We analyse the viability of this model by studying 
features of the background model and properties of the adiabatic UDM perturbations, 
which depend on the effective speed of sound and the functional form of the Jeans scale. 
As a result, theoretical constraints on the parameters of the models are found that allow 
for a behaviour compatible with observations. 
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1 Introduction 

The acceleration of the expansion of the Universe is, within the homogeneous and isotropic 
paradigm of cosmology, a well-accepted and observationally well-supported reality. The 
simplest possible framework explaining this acceleration is provided by the concordance 
ACDM model [1], where a cosmological constant A in Einstein equations sources the accel- 
eration, while Cold Dark Matter (CDM) is the main component for structure formation. 
Alternatives to A are various forms of Dark Energy (DE) [2] or a modified gravity theory 
[3]. A different approach is to consider models of Unified Dark Matter (UDM) where a 
single matter component is supposed to source the acceleration and structure formation at 
the same time (see [4, 5] and references therein). 

In a previous paper [4] the concept of UDM models with fast transition was introduced 
(see [6] for DE models with a sharp transition in their equation of state). In essence, these 
models are based on the idea that the Universe may have undergone a transition between 
a standard matter-like era, well described by an Einstein-de Sitter (EdS) model, and a 
ACDM-like epoch. If the transition is slow the differences with ACDM are negligible while 
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if the transition is fast, these models show interesting features 1 [4, 7]. The transition 
can be quantified as fast by looking both at parameters that govern the evolution of the 
background as well as at quantities that dictate the dynamics of perturbations. This 
analysis has been carried out in [4] for a specific barotropic model and will be generalised in 
this paper to a new class of models, but in essence the transition needs to be fast because: 
a) we are especially interested in background models that, at least in principle, can be 
clearly distinguished from a standard ACDM; b) as shown in [4], otherwise the evolution 
of perturbations is such that observational constraints are violated, in particular causing a 
strong deviation from the Integrated Sachs- Wolfe (ISW) effect occurring in ACDM models. 

The UDM models with fast-transitions may be an interesting alternative to other 
explanations of the observed acceleration of the Universe, providing a good fit to standard 
observables such as the CMB and the matter distribution [4] that require consideration 
of the perturbations, while at the same time showing interesting new features, leading to 
measurable predictions. So at least in principle, UDM models with a fast transition are 
able to avoid the fate of UDM models such as the generalised Chaplygin gas, which need to 
become indistinguishable from ACDM in order to survive observational tests, which spells 
their end [8], cf. [9, 10]. 

In contrast with more standard CDM+DE models, where the CDM component is 
perturbed and leads to structure formation while the DE component takes care of the ac- 
celeration of the background, with small or negligible effects on perturbations, the single 
UDM component must accelerate the Universe and provide acceptable perturbations. In 
particular, while CDM density perturbations evolve in a scale-independent fashion, this is 
not the case for UDM. In view of testing models against observations, e.g. with Markov 
Chain Monte Carlo methods and likelihood analysis, these differences may become com- 
putationally expensive. 

In view of this, and lacking a fundamental model, it is therefore essential to consider 
simple phenomenological models of the fast-transition paradigm that can be easily and 
efficiently implemented into numerical codes such as CAMB [11] and CLASS [12], and 
for which as much theoretical understanding as possible can be extracted from analytical 
calculations. Our first goal here is therefore to look at simple phenomenological recipes, 
such that the most important variables required for the numerical problem can be expressed 
analytically. We then introduce a specific simple model and analyse its properties in some 
detail, establishing what its range of validity should be if compared if observations. 

The rest of the paper can be outlined as follows. In section 2 we introduce the basic 
equations describing the background and the perturbative evolution for a general UDM 
model. In section 3 we explore three possible prescriptions for the dynamics of the UDM 
component. In section 4 we introduce a new UDM model with fast transition and study its 
background evolution, comparing it to a ACDM. In section 5 we analyse the properties of 

1 In a ACDM model of course there is an early matter era when A is negligible, with a smooth slow 
transition to the epoch when A plays a role. We are instead interested in models where there is a longer 
matter era that almost suddenly ends into a ACDM-like late behaviour; in other words, models where the 
UDM component evolves for long time in a CDM-like fashion, and suddenly a A-like term is switched on 
in the dynamics. 
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perturbations in this model, focusing on the evolution of the effective speed of sound and 
that of the Jeans length during the transition. The conclusions are drawn in section 6. 



2 Generalities of UDM models 

2.1 The background 

We assume a spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology, al- 
though many of the considerations in the following sections would apply in a more general 
context. The metric then is ds 2 = —dt 2 + a 2 (t)5ijdx' l dx^ , where t is the cosmic time, 
a(t) is the scale factor and 5{j is the Kronecker delta. The total stress-energy tensor is 
that of a perfect fluid with energy density p and pressure p, with its four-velocity: 
Tu,v = (p + p) u^Uy + pg^y Starting from these assumptions, and choosing units such that 
8ttG = c = 1 and signature { — ,+,+,-1-}, Einstein equations imply the Friedmann and 
Raychaudhuri equations: 

* -(!)'-§• <"> 

Oj 1 , S , v 

- = --p + 3p, 2.2 
a o 

where H = a/a is the Hubble expansion scalar and the dot denotes derivative with respect 
to the cosmic time. Assuming that the energy density of the radiation is negligible at the 
times of interest, and disregarding also the small baryonic component, p and p represent 
the energy density and the pressure of the UDM component. 

Independently from Einstein equations, projecting the conservation equations T^ v ]V = 
along one obtains the energy conservation equation 

p = -3H (p + p) = -3Hp (l + w) , (2.3) 

where w = p/p represents the equation of state (hereafter EoS) that is needed to close the 
system and is the quantity that characterises the background of our UDM model. 

When needed, we shall introduce different components, each with energy density pi. 
From this, we can define the dimensionless function 

0,(0) = , (2.4) 

Per [a) 

where p cr = 3H 2 is the critical density; values today will be denoted by the parameter O^o- 

2.2 Perturbations 

Assuming a perfect fluid, perturbations of the FLRW metric in the longitudinal gauge read 

ds 2 = -a 2 (r]) [(1 + 2$) dr] 2 - (1 - 2$) <%dxW] , (2.5) 

where $ represents the analogous of the Newtonian gravitational potential and we now are 
using conformal time rj. 
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Defining 

u = - 2.6 

and linearising the 0-0 and 0-i components of Einstein equations, one obtains the following 
second order differential equation for the Fourier component of u [4, 13-15]: 

d 2 u , 99 1 d 2 6 , _ N 

+ k 2 c 2 u - -—,u = , (2.7) 



drj 2 s 9d0 2 

where 



a +^), (2.* 



3(p + p) 

with z the redshift, 1 + z = a -1 . The quantity Cg in (2.7) is the effective speed of sound and 
characterises the perturbative dynamics of our UDM model, being also crucially involved 
in the growth of the overdensities 5p. Assuming adiabatic perturbations this is the same 
as the adiabatic speed of sound: 

drj 

Starting from Eq. (2.7), let us define the squared Jeans wave number: 

1 d 2 9 



c 2 e de 2 



(2.10) 



its reciprocal defines the squared Jeans length: Aj = a 2 /k 2 . 

An important aspect of UDM models is the possible manifestation of an effective 
sound speed significantly different from zero at late times: this generally corresponds to 
the appearance of a Jeans length (or sound horizon) below which the dark fluid does not 
cluster (e.g. see [4, 16, 17]). This causes a strong evolution in time of the gravitational 
potential, which at small scales starts to oscillate and decay, with effects on structure 
formation. In general, UDM models may also exhibit a strong ISW effect [15]. 

Thus, the squared Jeans wave number plays a crucial role in determining the viability of 
a UDM model, because of its effect on perturbations, which is then revealed in observables 
such as the CMB and matter power spectrum [4, 17]. As discussed in [4] there are two 
different regimes of evolution, respectively for scales much smaller and much larger than 
the Jeans length. In practice, any viable UDM model should satisfy the condition k 2 S> k 2 
for all the scales of cosmological interest, in turn giving an evolution for the gravitational 
potential <3? in Fourier's space of the following type (we are dealing with the gravitational 
potential after recombination so there is no more speed of sound due to radiation): 



H{rj) 
a(ri) 



a(fj) 2 dfj 



(2.11) 



The integration constant = <I> (0, k) T m (k) is fixed during inflation by the primordial 
potential $ (0, k) at large scales; T m (k) is the matter transfer function, describing the 
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evolution of perturbations through the epochs of horizon crossing and radiation-matter 
transition, see e.g. [18]. 

The explicit form of the Jeans wave number is 



2 pa 



,(l + w) 



dc, 



w 



s , 3(q 



w 



2(q 



w 



1 

+ 3 



(2.12) 



dp 6(1 + w) 

It is therefore clear from this expression that, if we want an analytic expression for kj in 
order to obtain some insight on the behaviour of perturbations in a given UDM model, we 
need to be able to obtain analytic expressions for p, p, w and c^. Unfortunately, it is not 
possible to find such expressions as functions of r/ (or t), simply because this requires the 
knowledge of an analytic expression for the scale factor as function of time, i.e. to solve 
the Friedmann equation (2.1), which in general is only possible for very special cases, as is 
well known. 



3 Three possible prescriptions 

A way out of this problem is to disentangle the evolution of the quantities of interest 
(p, p, w and (? s ) from Einstein equations, noticing that we can obtain these quantities as 
functions of the scale factor a if we use only the conservation equation (2.3). This has also 
the advantage that the expressions so obtained will be the same in any theory of gravity 
that satisfies the conservation equations. Eq. (2.3) becomes, using a as time variable: 

p' = --(p + p), (3.1) 

a 

where a prime indicates derivative with respect to a. 

We now briefly look at three different possible ways to prescribe the dynamics of the 
UDM component and to derive analytic expressions for the needed variables. 

3.1 Starting from w(a) 

Suppose that p/p = w is pre-assigned as a function of the scale factor: w = w(a). For 
instance, w(a) = wq + w a (l — a) is a typical phenomenological assumption well- motivated 
when setting observational constraints on dark energy models [19, 20]. In principle this 
is a convenient practical prescription to model UDM, because we have a good idea about 
what type of w(a) we should have. Then the adiabatic speed of sound 



s dp & of 

da 



can be computed from the conservation equation (3.1). Indeed, from p = wp we can 
compute p', then substituting the latter and p' from (3.1) to obtain 



aw' 



c 2 „ = w ; r . (3.3) 

3(1 + to) V ; 

This prescription however doesn't lead to analytic expressions for p(a) and p{a) in general, 
unless J 1+ ^( a ) da is integrable. 
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3.2 Starting from p(a) 

Prescribing the pressure as a function of the scale factor, p = p(a), can be useful, e.g. if one 
is dealing with a scalar field, in which case this is equivalent to prescribing a Lagrangian, 
see e.g. [5, 21, 22] and Refs. therein. Again, we have a good idea about the functional 
form that p(a) should have, so that this also seems a good starting point. Let's rewrite the 
energy conservation (3.1) as 

p' + -p = --p(a). (3.4) 
a a 

The homogeneous solution is p m oc a -3 , i.e. standard matter (dust), and for a given p(a) 
an analytic expression for p(a) can be found if E = 3 f a 2 p(a)da = J pdV is integrable, 
giving 2 p = E/V . With this prescription c 2 is immediately found, given p{a) and (3.4), 
but an analytic expression for w(a) can only be found if that for p(a) is found. 

3.3 Starting from p(a) 

It is perhaps less obvious what the functional form for p{a) should be, but some guess can 
be made in view of constructing UDM models with fast transition. In this case we want 
to recover a CDM-like behaviour, i.e. an EdS model, at early times (before the transition), 
i.e. p ~ p m = PM : oa~ s , and a CDM+DE behaviour after the transition. If we want to 
recover the simplest case of a ACDM at late times, we would have p ~ p\ + pm,oo -3 . 

Assuming that p(a) is given, we can compute p'(a) and then use (3.4) to obtain p(a) 
and then p'(a). In this case therefore any choice of p(a) guaranties analytic expressions for 
the EoS w(a) and the adiabatic speed of sound c 2 (a). 

To summarise, given a function (at least of class C 3 ) p = p(a) for the energy density, 
we have the following expressions for the quantities that enter into the Jeans wave number 
(2.12): 

a p' , 
w = --^-l, 3.5 
3 p 

2 a P" 4 ,„ (,\ 

C * = -37~3' (3 - 6) 



2 



dc 



dp 3/3' 



112 

III i // P 

ap + p - a— - 
P 



(3.7) 



Substituting in Eq. (2.12) from (3.5)-(3.7) we can obtain an analytic expression for the 
function A; 2 (a). Armed with this, we can obtain some insight about the behaviour of 
adiabatic perturbations in a model with a specified energy density p = p(a). 

4 Phenomenological UDM models with fast transition 
4.1 An overidealised model 

Before looking at a possible model for a UDM with fast transition, we now introduce an 
overidealised model, using a Heaviside function to describe an instantaneous transition. 



2 The expansion is adiabatic, and p m oc V \ thus we can interpret E as the energy of the system in the 
volume V . 
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This model can't serve our purposes, because it is clear from the expressions above that we 
need the function p(a) to be at least of class C 3 , but it is useful to get an idea of what we 
want to obtain. We assume that the Universe is well described by an EdS model before the 
transition, while for generality we describe the post-transition era with an "affine" model 
[17, 21, 23-25]: 



P 



(4.1) 



Pt ) a < a t 

.PA + (A"PA)(f) 3(1+a) a>a t 
Of course, any other post-transition model could be chosen. Here pt is the energy scale 
at the transition and p\ is the effective cosmological constant: 1 + a > and the energy 
density at late times tends to p\, so that the late time evolution in these models is a-la 
de Sitter, even if there is no cosmological constant in the Friedmann and Raychaudhuri 
equations, see [17, 23-25]. On the other hand, at is not an independent parameter: using 
the Friedmann equation (2.1) and the second of Eq. (4.1), and neglecting radiation, we get 



i - n A , 



n 



A,0 



3(l+c0 



(4.2) 



1. We can then 



where f^o/^A.o = Pt/ PA] this gives a redshift for the transition zt = a t 
interpret our UDM model after the transition as made up of the effective cosmological 



Pt ~ PA," 



constant p\ and an evolving part, with energy density at the transition p n 
that decreases after the transition. 

In the affine model the energy density of UDM is given by the second of Eq. (4.1). If 
one assumes this model all the way deep into the radiation era, the evolving part must have 
Wm and (? s — > in the past in order to recover standard matter domination at early 
times, thus the evolving part can be interpreted as the "dark matter" component today. 
Indeed, the value of the parameter a is extremely constrained if one assumes an affine model 
with adiabatic perturbations all the way back to recombination and beyond, a ~ 10 [17] 
(cf. also [26]), making this model indistinguishable from ACDM. 4 However, such a strong 
bound mainly comes from the matter power spectrum at small scales, k > kj, where the 
growth of perturbations is affected by a non- vanishing kj. Since this is an integrated effect, 
it is reasonable to expect that, if the affine-like evolution only starts below a transition red- 
shift zt, the bound will be much weaker. 

It is useful to explicitly incorporate a Heaviside function H(a — at) in Eq. (4.1), so 
that: 

'a t \3(i 



Pt 



f a t\ 3 , \ /' 

[—) + PA + (Pt — Pa) y- 



Pt 



at\ 3 
a / 



H(a - a t ) 



(4.3) 



For a = this reduces to 



P 



Pt[~) +PA 



1- - 



OtV 



~a) 



H(a - at) 



(4.4) 



3 In this two component interpretation the two energy densities satisfy their own conservation equations 
and therefore there is no coupling between them. 

4 The ACDM can be obtained as a sub-case of the UDM affine model for a = [23-25]; the constraint 
is weaker if the perturbations are not adiabatic, see [17]. 
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Figure 1. Illustrative parametric plot of p/ pa as a function of p/pa with pt/pA = 100 (z t ~ 5.34) 
for the model specified by Eq. (4.5). The ACDM model is represented here by the solid horizontal 
line p/pa = — lj while the line p = represents an EdS model, i.e. pure CDM. All models 
asymptotically evolve towards the effective cosmological constant pa- 

representing a sudden transition to ACDM. In the following, we shall restrict our attention 
to this sub-class of models. 

It is now clear that, simply replacing H(a — at) with a smoother transition function 
Ht(a — at), we can obtain simple UDM models with a fast transition. 

4.2 A simple model for the background 

In this paper, among the many known continuous approximations to the Heaviside function 
[27], we shall consider the only one that we found compatible with having (? s > 0: 



where the parameter ft represents the rapidity of the transition. In addition to at and ft, 
there is a third parameter in the model, which is p\, or, equivalently, the corresponding 
density parameter Oao- From inserting Eq. (4.5) in Eq. (4.4) we see that asymptotically 
in time, in the limit a — > oo, p — > p\, which implies p — > —p\. As already mentioned 
above, p\ plays the role of an effective cosmological constant, i.e. it is an attractor for 
Eq. (3.1). The Universe necessarily evolves toward an asymptotic de Sitter phase, i.e. a 
sort of cosmic no-hair theorem holds (see [28, 29] and refs. therein and [23-25]). In Fig. 
1 we show a parametric plot of p vs p, normalised to p\, where we have assumed, as an 
example, that the transition takes place at zt = 5.34, for a representative choice of values 
of ft. It can be seen that all models gradually approach the effective cosmological constant 




(4.5) 



PA- 
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In order to make our UDM model close to ACDM at late times, the transition must 
occur at relatively high redshifts, such that p t is quite larger than p\, which corresponds 
to a minimum value of the redshift zt- Otherwise, it could be difficult to have a good fit 
of supernovae and ISW effect data [4]. For instance, we need zt > 2.65 if we want to have 
pt > 20/3A- In Fig. 2 the evolution of the EoS w is depicted as a function of the scale factor 
for different values of pt/pA and /3. As shown, models with a smaller rapidity (3 have a 
background evolution more similar to that of the ACDM model at all times. On the other 
hand, a larger (3 implies a sharper transition between the CDM-like phase and the ACDM 
phase. Likewise, it clearly illustrates that the transition has to take place far enough in 
the past, i.e. pt is larger than p\, in order for the late time evolution of w to be close to 
that of the ACDM model. 

Finally, let us consider the evolution of the three components p r (radiation) p m and 
p\, represented by the dimensionless functions ili(a), Eq. (2.4). The evolution of these 
functions is shown in Fig. 3 for UDM models with fast transition (/5 = 500) occurring 
at different times, contrasted with pure the CDM (EdS) and ACDM models. For UDM 
and ACDM, it is assumed that today Oa.o = 0.72. In a flat Universe, £j J7j = 1, and 
when a component j dominates O,- ~ 1, with the other Oj ~ 0. It can be seen from 
these figures that if the transition is too late (today in the extreme case zt = 0) then 
the matter-radiation equality of the UDM model is basically the same as in a pure CDM 
model, i.e. much earlier than in ACDM. In addition, the effective cosmological constant of 
the UDM becomes dominant at a later time than in ACDM. Since the matter-transition 
equality dictates when matter perturbations inside the horizon start to grow, while the 
late dominance of the cosmological constant slows down this growth, the matter power 
spectrum in fast transition UDM models with too late a transition will be at odds with the 
observed one, with too much power on small scales. Conversely, we see from Fig. 3 that 
if the transition is at z% ~ 1 the matter-radiation equality is closer to that of ACDM, and 
it essentially coincides with the latter for an even earlier transition, Zt ~ 2. Basically, as 
long as the transition is at a redshift zt higher than the one at which f^A = in ACDM, 
the late time evolution of £Ia an d £l m in the UDM models is the same as in ACDM. 

4.3 Angular diameter distance 

The angular diameter distance is an important quantity that comes into play in current 
observations of weak lensing, BAO and galaxy clustering and will become even more im- 
portant for comparing models against the new data that will become available from surveys 
such as DES, Planck and Euclid. It is also relevant to measurements of CMB anisotropics. 
In view of this, we now briefly comment on the deviation of the angular diameter distance 
in our UDM model from that of ACDM. The angular diameter distance, d,A, is defined 
as the ratio of the actual size Ax of an object and the angle A6 this object subtends 
orthogonal to the line of sight and can also be expressed as: 




(4.6) 
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where H{z) is the Hubble function. In Fig. 4, we have show the angular diameter distance 
normalised with the Hubble distance, dn = H^ 1 , for the ACDM model and our UDM 
model with a fast transition with f3 = 500. As it can be seen, the larger the transition 
redshift, the smaller the departure of dA for the UDM with respect to the ACDM is. Notice 
that dA does not increase indefinitely as z — > oo; it turns over at z ~ 1.5 and thereafter 
more distant objects actually appear larger in angular size. In Fig. 5 the difference between 
dAjdu for the ACDM and for the UDM with a fast transition at /3 = 500 is depicted. 



5 The Jeans scale and the gravitational potential 
5.1 The Jeans wave number 

We now focus on the Jeans wave number for our UDM model and investigate its behaviour 
as a function of the speed of sound, in particular around p = pt, which corresponds to the 
middle of the transition where the speed of sound is at its peak. 

By inspection of Eq. (2.12) we see that a large k 2 can be obtained not only when 
c 2 — > 0, but when c 2 changes rapidly as well. In other words, when Eq. (2.12) is dominated 
by the p dc 2 /dp term we may say that the EoS is characterised by a fast transition. These 
two quantities, c 2 and p dc 2 /dp, are depicted in Fig. 6 as functions of p/pa- 

Thus, viable adiabatic UDM models can be constructed which do not require c 2 <C 1 
at all times if the speed of sound goes through a rapid change, a fast transition period 
during which k 2 can remain large, in the sense that k 2 -C k 2 for all scales of cosmological 
interest to which the linear perturbation theory of Eq. (2.7) applies. 

When we consider a fast transition, it is interesting to compare the term pdc 2 /dp with 
the remaining ones contained in the squared brackets of Eq. (2.12) for the Jeans wave 
number, that is: 

R 1 f2 , , 3(c s 2 - w) 2 - 2{cl - w) 1 

In Fig. 7 we plot pdc 2 /dp, B and the Jeans wave number kj as functions of p/p\- 
We have chosen p t = 100pA i n order to consider a transition sufficiently back in the Dark 
Matter epoch and a value of (3 = 500 in order to show a fast transition. 

From Fig. 7 we learn that p dc 2 / dp is negative for p > p t , then for p ~ p t it increases 
becoming positive and intersecting the B curve a first time for p < p t . For smaller values of 
the energy density, pdc 2 /dp decreases again to zero, again intersecting the B curve. Since 
the difference between the two curves is very large, we have depicted them in the right top 
panel of Fig. 7 choosing a logarithmic scale. For this reason, the negative part of pdc 2 /dp 
has been omitted. It can be also seen that the relative maximum of the Jeans wave number 
between the two zeros of the curve approximately corresponds to the point where dc 2 /dp 
reaches its maximum value. 

The place where the curves p dc 2 /dp and B intersect correspond to the vanishing points 
of the Jeans wave number fcj, as it can be seen in the bottom panel of Fig. 7. In general, 
around these points the corresponding Jeans length becomes very large, possibly causing all 
sort of problems to perturbations, with effects on structure formation in the UDM model. 
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Defining c 2 as in Eq. (3.6), in Fig. 8 we plot the Jeans wave number as a function of the 
redshift for various values of j3 and pt/ph- In this figure, for sufficiently high (3 we note 
that i) in general the Jeans wave number becomes larger, with a vanishingly small Jeans 
length before and after the transition, and ii) it becomes vanishingly small for extremely 
short times, so that the effects caused by its vanishing are sufficiently negligible, as we are 
going to show in the next subsection when we analyse the gravitational potential <5. 

To conclude, we shall make some comments on building phenomenological UDM (or 
DE) fluid models intended to represent the homogeneous FLRW background and its linear 
perturbations. A fast transition in a fluid model could be characterised by a value of c 2 » 1 
during the transition. However, as the FLRW background is homogeneous, c 2 = p/p does 
not actually represent a speed of sound because there is nothing that could propagate in 
this case. Regarding the linear perturbations, for the scales in which k <C kj the solution 
of the Eq. (2.7) for the gravitational potential is Eq. (2.11) and therefore for such scales 
there is no superluminal propagation. The reason behind this is that Eq. (2.7) is the 
Fourier component of a wave equation with potential 6"/6, and the latter does not allow 
propagation for k <C kj. So, we can always build the fluid model in such a way that all 
scales smaller than the Jeans length A <C Aj correspond to those in the non-linear regime, 
for which this model may not apply. In order to study the behaviour of the perturbations 
of a UDM model at these scales, we would have to go beyond the perturbative regime 
investigated here. That would possibly imply to build a more refined fluid model that could 
maintain causality and at the same time be able to deal with the increased complexity of 
small scale non- linear physics. 

5.2 The gravitational potential 

The differential equation that governs the behaviour of the gravitational potential <£> in our 
model in terms of the scale factor a is obtained from Eq. (2.7): 



d 2 $(k,a 



da 2 



,H da a a J da 



2 dU 1 



.21,2 



cik 



+ a2 ( 1 + 3c?) + a 



$(k,a) = , 
(5.2) 

where T~L = ^/a is the conformal time Hubble function. Also, T~L = aH, c 2 is given in 
Eq. (3.6) and we have assumed plane-wave perturbations <3?(x, a) oc <5(k, a) exp (ik • x) of 
comoving wave- number k = |k|. 

On the other hand, the ACDM gravitational potential, <1?a> solves Eq. (5.2) for c s = 0: 



d 2 <S> 



K / 1 dH | 4\ d$ A | 



2 dU 1 

aH da a 2 



$A = 



(5.3) 



da 2 da ' a J da 

where 

V 2 = Hi (O A , a 2 + n^oa- 1 ) , (5.4) 

if we assume that the energy density of radiation is negligible and ignore the contribution 
of the baryon matter. 

Comparing Eqs. (5.2) and (5.3) we see that the gravitational potential in our UDM 
model has the same evolution as in a ACDM universe in the limit when c s — > 0. In 
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particular, both models behave as an EdS model at early times, so that for UDM c s — > 
and in both models H 2 ~ a -1 , and $ is constant in this regime, as in an EdS universe, as 
is well known. 

The normalised initial conditions are $A(k; a rcc ) = 1 and d^\/da\ a ^ c = 0, where a rcc 
stands for the scale factor at recombination time. Since the class of UDM models we 
consider here is constructed to behave as the ACDM model in the early Universe, we thus 
set the same initial conditions for both the UDM and the ACDM gravitational potentials. 

In ACDM, the background evolution causes a gradual time evolution of the gravita- 
tional potential when the cosmological constant starts to dominate [30]; this causes an ISW 
effect. On the other hand, in our UDM model the evolution of the gravitational potential is 
determined by the background and the perturbative evolution of the single dark fluid and, 
crucially, by the adiabatic speed of sound c 2 s . The gravitational potential stays constant 
before the transition, during which a sudden rapid evolution of <1> is induced. The subse- 
quent evolution is in general scale dependent: for scales k > kj, <E> oscillates and decays; for 
larger scales, k < kj, the evolution of $ becomes scale independent and is governed by the 
evolution of the background, mainly by H and to a small extent by c^, and <3? approaches 
its ACDM behaviour in a way that depends mainly on the rapidity /3 and from the epoch 
of the transition, zt- In general, we expect an ISW effect starting from the transition which 
can be very different from the ACDM one, cf.[4]. We show the behaviour of $ in Fig. 9, 
where we explore its dependence on the background parameters /3 and z t (or, equivalently, 
at)- 

We already know from the evolution of w in Fig. 2 that for (3 < 200 the transition 
is not that fast. This is also apparent in Fig. 9, where we have plotted the normalised 
gravitational potential &k{ z ) = $(k; z)/<]?(0; 10 3 ) as a function of the redshift z for k = 0.2 
h Mpc -1 and different values of (3 and z t . For wavenumbers k > k n i — 0.2 h Mpc -1 , we 
expect non-linear matter over-densities contributions to the evolution of the gravitational 
potential to become important. 

As expected, for large enough values of /3 the gravitational potential is practically 
constant in time, corresponding to a pure matter (EdS) background evolution, until z ~ 
zt. From z < z t onwards, its value decreases very fast until it finally approaches the 
gravitational potential of the ACDM model. The expected strong ISW effect caused by 
this behaviour would be mostly due to the particular expansion history and can be very 
different from the ACDM one. As a matter of fact, these differences are smaller for z% > 2 
and become smaller and smaller at higher transition redshifts, until for zt ~ 100 and larger 
the fast transition UDM models practically become indistinguishable from ACDM cf. [4]. 

As we have already pointed out, our UDM model allows the value w = — 1 for a — > oo, 
i.e, it admits an effective cosmological constant energy density p\ at late times. Hence, 
if we wanted to compare the predictions of our UDM model with observational data, we 
would follow the prescription given in [4], where the density contrast is 5 = 5p/pA and 
PA — p — PA is the clustering "aether" part of the UDM component [23, 31]. In UDM 
models gravity is described by General Relativity, but to link the density contrast with the 
gravitational potential at scales much smaller than the cosmological horizon we only need 
the Newtonian Poisson equation. For z < z rec , where z Tec is the recombination redshift 
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(z Tec ~ 10 3 ) we then have: 

PA 

To conclude, we have argued that for an early enough fast transition with f3 > 500 
and zt > 2 our UDM model should be compatible with observations. On the other hand, 
a study of the matter and CMB power spectra is needed to study the viability of models 
with 10 < /3 < 500, and those with /3 > 500 and z t < 2. We shall undertake this work in 
the future. 



6 Conclusions 

UDM models, when compared with the standard DM + DE scenario, specifically ACDM, 
are in principle interesting because the dynamics of the Universe can be described with a 
single component in the matter sector which triggers the accelerated expansion at late times 
and is also able to cluster and produce a satisfactory structure formation. The challenge 
for UDM models is however to satisfy observational constraints while maintaining features 
that can make then distinguishable from ACDM, otherwise they loose interest [8]. 

After considering (in Sec. 3) three possible prescriptions for building simple UDM 
models, in this paper we have introduced and examined a new class of UDM models with 
a fast transition between an early matter era and a later ACDM-like phase, building on 
previous work [4, 7]. In particular, here we have focused on models in which all the variables 
of interest can be expressed analytically, so that in principle they could be efficiently 
implemented into numerical codes such as CAMB [11] and CLASS [12]. 

A general feature of UDM models is in the possible difference of the expansion history 
with that of ACDM, causing, among other features, a strong ISW effect incompatible 
with observations. In addition, in UDM models the effective speed of sound may become 
significantly different from zero. This corresponds, in general, to the appearance of a Jeans 
length (or sound horizon) below which the dark fluid cannot cluster and which, if large 
enough, can cause a strong evolution in time of the gravitational potential, preventing 
structure formation at small scales. In building satisfactory UDM parametric models it 
is therefore crucial to find the region in parameter space where the Jeans length remains 
small enough, well beyond the linear regime that we explore here. 

UDM models with a fast transition, first introduced in [4, 7] are a viable and interesting 
alternative to ACDM because they seem to survive observational tests while maintaining 
interesting features. The general phenomenological UDM models we have introduced in 
Sec. 4 are characterised by a fast transition between a standard matter era and a post- 
transition epoch described by an "affine model" [17, 21, 25] with affine parameter a. We 
have then focused on the a = case, which represents a sudden transition to ACDM. 
In this section we have also compared the angular diameter distance between ACDM and 
our UDM with fast transition, finding small differences of the order of percent when the 
transition is fast enough. 
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Finally, in Sect. 5, in order to study the viability of our UDM model, we have carried 
out a study of the functional form of the Jeans scale in adiabatic UDM perturbations. 
In doing so, we have found analytical expressions for the quantities involved in the Jeans 
wave number and have shown that our model presents a small Jeans length even when a 
non- negligible sound speed is present. Subsequently, we have analysed the properties of 
perturbations in our model, focusing on the evolution of the effective speed of sound, the 
Jeans scale and the gravitational potential. In general, in building a phenomenological 
model, we have chosen its parameter values in order to always satisfy the condition fccfcj 
for all k of cosmological interests to which linear theory applies. In this way, we have been 
able to set theoretical constraints on the parameters of the model, predicting sufficient 
conditions for the model to be viable. We have argued that for large enough values of the 
rapidity [5 of the transition and z% our model should be compatible with observations. 

Computing the CMB and the matter power spectra for our model, for a wide range of 
parameters values, as well as a full likelihood analysis for this model and its parameters, 
including will be the subject of a forthcoming work. 
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Figure 2. Evolution of the UDM parameter w = p/p as a function of the scale factor a for 
Pt/PK — 10 (a t ~ 0.35, z t ~ 1.85) (top panel) and pt/pA = 20 (a t ~ 0.27, z t ~ 2.65) (bottom 
panel). For reference we also plot: the w = line representing a flat pure CDM model (an EdS 
Universe); the w — — 1/3 line representing the boundary between the decelerated and the accelerated 
phases and the thin black curved line representing the evolution of the total w for the ACDM model 
with f^o = 0.72. The higher pt/pA, the earlier the UDM w transits to that of a ACDM. In the 
future, for a > 1, all models evolve to a dc Sitter phase with w = — 1. 
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Figure 3. Density dependence of the various components for UDM models with a fast transition 
(j3 = 500) occurring at different times, contrasted with the pure CDM (EdS) and ACDM models, 
assuming f2 A ,o — 0.72 today. Top panel: z t = 0; medium panel: z t = 1; bottom panel: z t = 2. 
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Figure 4. The dimensionless angular diameter distance (IaI&h for the ACDM and the UDM model 
with a fast transition at (3 = 500 for different transition redshifts. 



0.00 



-0.01 



dA (udm) 



d 



A (ACDM) 



1 -0.02 



-0.03 



-0.04 















\ 

s 


Zt 


= 5.25 




z t 


= 4 ; 




Zt 


= 2.7 : 


: yS = 500 











2 3 
z 



Figure 5. Difference between (Ia/cLm for the ACDM and for the UDM model with a fast transition 
at j8 = 500. 
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Figure 7. Left-top panel: evolution of p dc^/dp (solid line) and B (dashed line) as functions of 
p/pA- Right-top panel: evolution of p dc^/dp (solid line) and B (dashed line) as functions of p/pa 
in a logarithmic scale. Bottom panel: evolution of the Jeans wave number kj as a function of p/ pa 
where kj is in h Mpc -1 units. The choice of parameters is pt/pA = 100 (zt ~ 5.34) and (3 — 500. 
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Figure 8. The Jeans wave number kj (in h Mpc -1 units) for /3 = 10 (left-top panel), j3 = 50 
(right-top panel), /3 = 100 (left-bottom panel) and /3 = 500 (right-bottom panel). In each panel, 
the three dotted, dash-dotted and dashed lines correspond to pt/pA = 10,20,50, respectively; the 
black solid line corresponds to Pt/pA = 100. 
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Figure 9. Illustrative plots of the gravitational potential $(k; z) as a function of the redshift z for 
ACDM and for our UDM model for k = 0.2 h Mpc^ 1 and different values of j3 and z t . The black 
solid line corresponds to the gravitational potential in the ACDM model with £Ia,o = 0.72. 
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